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Abstract 

As an alternative to dark energy it has been suggested that we may be at the center of an in- 
homogeneous isotropic universe described by a Lemaitre-Tolman-Bondi (LTB) solution of Einstein's 
field equations. In order to test such an hypothesis we calculate the low redshift expansion of the 
luminosity distance Dl(z) and the redshift spherical shell mass density mn(z) for a central observer 
in a LTB space without cosmological constant and show how they cannot fit the observations implied 
by a ACDM model if the conditions to avoid a weak central singularity are imposed, i.e. if the matter 
distribution is smooth everywhere. Our conclusions are valid for any value of the cosmological con- 
stant, not only for f^A > 1/3 as implied by previous proofs that qQ PP has to be positive in a smooth 
LTB space, based on considering only the luminosity distance. 

The observational signatures of smooth LTB matter dominated models are fundamentally different 
from the ones of ACDM models not only because it is not possible to reproduce a negative apparent 
central deceleration (?g PP , but because of deeper differences in their space-time geometry which make 
impossible the inversion problem when more than one observable is considered, and emerge at any 
redshift, not only for z = 0. 



1 



I. INTRODUCTION 
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and the WMAP measurement 



High redshift luminosity distance measurements 
43] of cosmic microwave background (CMB) interpreted in the context of standard FLRW 
cosmological models have strongly disfavored a matter dominated universe, and strongly sup- 
ported a dominant dark energy component, giving rise to a pos itive cosmological acceleration. 
As an alternative to dark energy, it has been proposed 13|, |l4| that we may be at the center 
of an inhomogeneous isotropic universe described by a Lemaitre-Tolman-Bondi (LTB) solution 
of Einstein's field equations, where spatial averaging over one expanding and one contracting 
region is producing a positive averaged acceleration ap. Another more general approach to 
map luminosity distance as a function of redshift Dl(z) to LTB models has been recently pro- 



posed 



30 



3jJ ], showing that an inversion method can be applied successfully to reproduce the 



observed Dl(z). 

The main point is that the luminosity distance is in general sensitive to the geometry of 
the space through which photons are propagating along null geodesies, and therefore arranging 
appropriately the geometry of a given cosmological model it is possible to reproduce a given 
Dl(z). For FLRW models this corresponds to the determination of and fl m and for LTB 
models it allows to determine the functions E(r), M(r),tb(r). 

Another observable which could be used to constrain LTB models is the redshift spherical 
shell mass mn(z) {2], recently calculated Q] for a central observer up to the fifth order, which 
is the product of the number of sources n(z) times their mass m(z). This observable can be 
related to observations by using the following relation with the total redshift rest mass Mt{z) 
, which is defined as the total rest mass of all the sources with redshift equal or less than z : 



M T (z) = / ATimn(z')dz\ {I) 





A7tmn{z) = !^I1±L_ (2) 

dz 

Once an actual accurate estimation of mn(z) based on observational data from galaxy sur- 
veys will be available, this quantity could be used to distinguish without ambiguity between 
LTB and ACDM models, but an even better candidate as a discriminating observable would 



451 ]. which has the additional advantage 



be the redshift spherical energy E R ss introduced in 
of avoiding source evolution effects. 

Erss is in fact obtained by integrating mn(z) over varying redshift intervals AZ(z) 



z+AZ(z) 

E RSS {z)= I ^mn{z')dz' (3) 



t(z) -t(z + AZ{z)) = At (4) 

which preserve the constant time interval At. By choosing At to be sufficiently smaller than 
the time scale over which astrophysical evolution is important, we can avoid the effects of the 
source evolution. 

In this paper calculate the low redshift expansion of mn(z) and Dl(z) for flat ACDM and 
matter dominated LTB. We then show how, if the conditions to avoid a central weak singularity 
are imposed, it is impossible to mimick dark energy with a LTB model without cosmological 
constant for both these observables, giving a general proof of the impossibility to give a local 
solution of the inversion problem for a smooth LTB model. This central singularity is rather 
mild, and is associated to linear terms in the energy density which lead to a divergence of the 
second derivative, so non smooth LTB models could still be viable cosmological models. It can 



be shown that the inversion problem 34j can be solved if the smoothness conditions we are 



imposing are relaxed. This implies that the numerical solutions of the inversion problem which 
have been recently proposed [2^ must contain such a weak central singularity. 

Our proof of the impossibility of the local solution of the inversion problem is general, since 
we don't use any special ansatz for the functions defining LTB models, and it does not depend 
on the particular value of Q\. In this regard we obtain a much stronger result than the one 
already well known fact that qQ PP cannot be negative for a smooth LTB model. In fact even if 
qQ PP > 0, which for a flat ACDM would imply < ~, it would still be impossible to mimick 
the cosmological constant for both mn(z) and Dl{z). 

This means that the difference between a smooth spherically symmetric pressureless matter 
dominated Universe and a homogeneous Universe with a cosmological constant goes beyond 
the sign of qo PP , and corresponds to distinctive observational features not only at z — 0. From 
a geometrical point of view we can interpret our results as an example in which redshift space 
observations allow to distinguish between different cosmological models despite the fact that 
they don't allow us to directly probe the actual local geometry of the Universe but only to 
access some "compressed" information about it. 
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II. LEMAITRE-TOLMAN-BONDI (LTB) SOLUTION 



Lemaitre-Tolman-Bondi solution can be written as 



0Jio| 



ds 2 



-dt 2 + 



(R, r ) dr 



+ R 2 dQ 2 



(5) 



1 + 2E 

where R is a function of the time coordinate t and the radial coordinate r, R = R(t,r), E is 
an arbitrary function of r, E = E(r) and R, r = dR/dr. 
Einstein's equations give 



p(t,r) 



2E(r) + 2M(r) 



R 2 
2M, r 
R?R n 



R 3 



(6) 
(7) 



with M = M(r) being an arbitrary function of r and the dot denoting the partial derivative 
with respect to t, R = dR(t,r)/dt. The solution of Eq. (J6]) can be expressed parametrically in 
terms of a time variable r = /' dt'/R(t', r) as 

M(r) 



Y r,r) 



t{r,r) 



-2E(r) 

M(r) 
-2E(r) 



1 - cos ( yJ-2E(r)T 
1 



sin 



-2E(r) 



-2E(r)T 



+ t b (r) 



(9) 



where Y has been introduced to make clear the distinction between the two functions R(t, r) 
and Y(r,r) which are trivially related by 



R(t(T,r)) = Y(r,r) 



(10) 



and tb{r) is another arbitrary function of r, called the bang function, which corresponds to the 
fact that big-bang/ crunches can happen at different times. This inhomogeneity of the location 
of the singularities is one of the origins of the possible causal separation 32| between the central 
observer and the spatially averaged region for models with positive ao- 
We introduce the variables 

R(t,r) 2E(r) 6M(r) 



A(t,r) 



k(r) 



Po{r) 



(11) 



so that Eq. (jSJ) and the Einstein equations (JS} and (|7j) are written in a form similar to those 
for FLRW models, 



ds * = - dt * + A 1 



l + 



A, r T 



dr 2 



1 — k{r)r'' 



+ r 2 d£l\ 



121 



4 



(J) 

p{t,r) 



K r ) , Po( r ) 
' A 2 3A 3 

(Po^ 3 ),r 



3A 2 r 2 (Ar) 

The solution of Eqs. OH]) and (Q can now be written as 

Po(r] 



(13) 
(14) 



afar) 
tfa r ) 



6k(r) 

Po(r) 
6k(r) 



cos 



7] 



k{r) 7] 



sm 



k(r) 



k(r) i] 



+ t b (r) 



(15) 
(16) 



Furthermore, without loss of 
Po = constant, corresponding 



where r/ = rr = f* dt' / 'A(t',r) and A(t(r), r), r) — a(r), r). 

In the rest of paper we will use this last set of equations 
generality, we may set the function po{r) to be a constant, po(r] 
to the choice of coordinates in which M(r) oc r 3 , and we will call this, following [13|. the FLRW 
gauge. 

We need three functions to define a LTB solution, but because of the invariance under general 
coordinate transformations, only two of them are really independent. This implies that two 
observables are in principle sufficient to solve the inversion problem of mapping observations 
to a specific LTB model, for example the luminosity distance Dl(z) and the redshift spherical 
shell mass m(z)n(z) = mn(z). 

As observed by 2£j] , there has been sometime some confusion about the general type of 
LTB models which could be used to explain cosmological observations, so it is important to 
stress that without restricting the attention on models with homogeneous big bang, a void is not 
necessary to explain both Dl(z) and m(z)n(z) = mn(z) without the cosmological constant, but 
we will prove that this is only possible for LTB models that contain a weak central singularity. 



III. GEODESIC EQUATIONS 



We will adopt the same method developed in 17| to find the null geodesic equation in the 
coordinates {i],t), but here instead of integrating numerically the differential equations we will 
find a local expansion of the solution around z = corresponding to the point (to, 0) = fa, t), 
where t = tfa,r). We will also provide more details about the geodesic equation derivation 
which were presented in |l7| in a rather concise way. We will indeed slightly change notation to 
emphasize the fully analytical r.h.s. of the equations obtained in terms of (77, t), on the contrary 
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of previous versions of the light geodesic equations which require some numerical calculation of 
R(t,r) from the Einstein's equation ([6]). 

For this reason this formulation is particularly suitable for the derivation of analytical results. 

The luminosity distance for a central observer in a LTB space as a function of the redshift 
is expressed as 

D L (z) = (1 + zfR (t(z),r(z)) = (1 + zfr{z)a (v(z),r(z)) , (17) 

where ft(z),r(z)) or ((r)(z), r(zfj is the solution of the radial geodesic equation as a function 
of the redshift. 

The past-directed radial null geodesic is given by 

dT(r) -R tr (t,r) . . 

— — = /(T(r),r); f{t,r) = —= 18 

d r Jl + 2E r) 



where T(r) is the time coordinate along the null radial geodesic as a function of the the coor- 
dinate r. 

From the implicit solution, we can write 

T(r) = t(U{r),r), (19) 
dT(r) = MdMrf + dt 
dr dr) dr dr' 

where U(r) is the rj coordinate along the null radial geodesic as a function of the the coor- 
dinate r. 

Since it is easier to write down the geodesic equation in the coordinate (t, r) we will start 



from there 



12|: 



dr y/l + 2E(r(z)) 

dz {l + z)R>[r(z),t{z)}' 
dt R'[r{z),t{r)] 
dz ~ (l + z)R'[r(z),t(z)]' 

(22) 

where the ' denotes the derivative respect to r and the dot ' the derivative respect to t. These 
equations are derived from the definition of redshift and by following the evolution of a short 
time interval along the null geodesic T(r). 

The problem is that there is no exact analytical solution for R(t,r), so the r.h.s. of this 
equations cannot be evaluated analytically but requires to find a numerical solution for R first 



[29] , and then to integrate numerically the differential equation, which is a quite inconvenient 
and difficult numerical procedure. 

Alternatively a local expansion for R(t,r) around (io>0) corresponding to the central ob- 
server, could be derived and used in eq.( )22|) . but being an expansion will loose accuracy as the 
redshift increases. 

For this reason it is useful for many numerical and analytical applications to write the 
geodesic equations for the coordinates (r],r), 

dr] d r t(r],r) — F(r],r) 



dz (1 + z)d v F(i],r) 
dr a {Vi r ) 



dz (1 + z)d v F(r), r) 



p(Vi r ), (23) 
qfar), (24) 



F(r],r) = 1 =[d r (a(r),r)r) + d v (a(r],r)r)d r r]] , (25) 

Jl — k(r)r 2 

where rj = U(r(z)) and F(rj,r) = f(t(r],r),r). It is important to observe that the functions 
p, q, F have an explicit analytical form which can be obtained from a(r), r) and t(r), r) as shown 
below. 

The derivation of the implicit solution a(t],r) is based on the use of the conformal time 
variable rj, which by construction satisfies the relation, 

drj{t, r) 



dt 

This means 



a' 1 . (26) 



t(r),r) = t b (r)+ / a(r] ,r)dr) , (27) 



o 



dt = a(r), r)dri + (j* + t' b (r)^ dr , (28) 

In order to use the analytical solution we need to find an analytical expression for F and F^. 
This can always be done by using 



3 sin (r)\Jk(r)j — r] (2 + cos {^\Jk{r) \ \jk~{r)\ 



I2k{r-f/ 2 I.*""" V v '" v 'J ''V ' " w V' /V '" v V v " v VJ +t ' b ^' ^ 
d r rj = — a(?7, r) _1 9 r t (31) 

In this way the coefficients of equations ( !23l) and ( )24l) are fully analytical, which is a significant 
improvement over previous approaches. 



IV. CALCULATING D L (z) AND mn(z) 



Expanding the r.h.s. of the geodesies equation we can easily integrate the corresponding 
polynomial q(z),p(z), to get r(z) and f](z). It can be easily shown that in order to obtain 
Dl(z) to the third order and mn(z) the fourth we need to expand r(z) to the third order and 
7)(z) to the second. 



In the FLRW gauge in order to to avoid a weak central singularity 15[ we need the following 
expansion for k(r) and tb(r): 

k(r) = k + k 2 r 2 (32) 
t b (r) = t b + t\r 2 (33) 

which are based on the fact that taking only even powers the functions are analytical every- 
where, including the center. 

Here we will not give the formulae in terms of rj and trigonometric functions, since they 
are rather complicated and not relevant to the scope of this paper, but rather introduce the 
following quantities: 



t m tan(^p) 2 p 

«o = a{T} , 0) = 7= (34) 

3Man(^|2o) 2 + 3£; 

3k 3 J 2 (tan(^p) 2 + l 
H = V " 2 L (35) 

tan (v^fo)3 po 

go 4(tan(^) 2 + l) (36) 



where we have used 



#o = ^4 (37) 
a(t ,0) 

a(t ,0)d(t ,0) 
a{t ,0) 2 

The derivative respect to t is denoted with a dot, and is calculated using the analytical 
solution a(r],r) and the derivative respect to rj : 

a = d t a = d^aa^ 1 (39) 

We may actually set ao = 1 by choosing an appropriate system of units, but we will leave it 
in order to clearly show the number of independent degrees of freedom of the problem and to 
emphasize the difference with AC DM models. 
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It is important to observe that in general, without imposing the smoothness conditions, we 
have eight independent parameters 



Po ,r ]0 ,t b o ,t b 1 ,t b 2 ,k ,k 1 ,k2 (40) 

The conformal time coordinate of the central observer t]q is not really independent, in the 
sense that it has to be consistent with the age of the universe, which, if we assume the inhomo- 
geneities to be only local, should be approximatively the same as the one estimated in AC DM 
models. 

We get six constraints from the expansion of mn{z) and Dl(z) respectively to fourth an 
third order, so in principle without imposing the smoothness conditions t\ = k\ = we should 
be able to solve the inversion problem of locally mapping a LTB model to any given ACDM 
model, even taking into account the fact that, as we will show later, mn(z) and D L (z) don't 
depend on t b Q . 

It is also clear that a simple preliminary argument based on counting the number of inde- 
pendent parameters imply that the inversion problem cannot be solved for a smooth model, 
i.e. if t\ — k\ — 0, since we have six constraints and only five truly independent parameters. 

We will show this more explicitly in the rest of the paper. 

After re-expressing the results in terms of H ,q , a>o we get 

T)(z) = r] + Vi z + V2Z 2 
r(z) = r\z + r 2 z 2 + r 3 z 3 



V2 



V2^T (a 4 # 4 {Aql - 3q + 1) - 2a 2 H$(l - 2q ) 2 t b 2 - 3k 2 ) + 2k 2 (q + 1) arctan ^2qJ^T 

2a^(2g - 1)5/2 



a H ^ ^ 

r 2 = -7T-TT ( 43 ) 



2a H 
1 



^3 



2a 5 Hl(2q - 1 



1 3 



(2q - 1) (a^l - 2<7o) 2 (<?o + l) " 2a^ 3 (l - 2g ) V2 - ^ 2 q + h) + 

(44) 



+6/c 2 y2go - lq arctan y 2g - 1 

It is important to observe that this formulae are valid, by analytical continuation, also for 
negative k , i.e. also for q < 1/2. 
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We can then calculate the luminosity distance 



Dr(z) 



D\ = 4- 



Dl 



1 + z) 2 r(z)a(rj(z), r(z)) = D\z + D 2 L z 2 + D 3 L z 3 
1 

H 

i - go 

2Ho 



2aiHl(2q - 1)5/2 



2q - 1 (a 4 # 4 (l - 2g ) 2 (g - l)?o - 2a 2 H 3 (4g 3 - 3g + l) t\ - 9/ 



-Qk 2 q (qo + 1) arctan W2g - 1 



From the definition of mn(z) and the equation for the energy density we can write 



4:7imn(z)dz = pd 3 V 



from which by using dr = [dr/dz)dz we get 



8nM' 



Jl — k(r)r- 



--dr 



(49) 



mn(z) 



2M'(r(z) 



dr(z) 



p r{z) 



dr(z) 



yjl - k(r(z))r(z) 2 dz yjl - k(r(z))r(z) 2 dz 
where in the last equation we have used the FLRW gauge condition M(r) 
allows to calculate mn(z) directly from r(z). 
We finally get: 



(50) 



Por 3 /Q, which 



mn{z) 
mri2 



mn 2 z 2 + mn 3 z 3 + mn a z A 

6g 

H 

12q (q + 1) 



H 
3g 



2atHl(2q - 1) 



20a 2 H 3 (l - 2qofqot b 2 + 10fc 2 (l - 5g )) + Q0k 2 J2q - lq 2 arctan J2q - 1 



(2q - l)(a^ 4 (l - 2g ) 2 (15g i + 14g + 13) + 



(51) 
(52) 

(53) 
(54) 



As it can be seen the effects of large-scale inhomogeneities show only from the fourth order 



and while also the fifth order has been calculated 
need it. 



4j , we don't report it here since we don't 



V. RELATION WITH APPARENT OBSERVABLES 

In this section we will briefly show the results obtained so far with apparent observables: 
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H app {z) 



d (D T (z) 



dz \ 1 + z 



(55) 



Qapp{z) = Tz (t+t) = H ~ 1{z) ™ 

aPP{z) = _ x _ pmi (57 ) 

H v ' dln{\ + z) v ' 



By apparent here we mean that they are obtained using the relations expressed in eq. ( 15511571 ) , 
which are true only in a flat FLRW space, but assuming that the observed luminosity distance 
used in these formulas correspond to a LTB model. 

As it could be expected on the basis of simple continuity arguments, for a smooth LTB 
models the central values of the apparent observables coincide with their dynamical definitions 
from eq.dSHMD : 

H = 4Mv| = H app (0) = fl?* (58) 
a(t ,0) 

g = - " (t °:° ) ^°; 0) = = m 

a(t ,oy 

While it is worth mentioning this relation in order to better understand the physical meaning 
of the quantities we have introduced, it is important to observe that eq. fl34H35|) can be just 
considered a convenient and natural re-definition of the parameters of the problem based, and 
that we don't't need neither use anywhere in our calculations the concept of apparent q app (z) 
and H app (z) introduced above. 

It is also important to observe that our results are consistent with the fact that qQ PP must be 
positive in absence of a weak central singularity, since the value calculated in eq. (|36l) is always 
positive for any value of 7] or k . It can be shown that equations (I58H59P are not satisfied for 



not smooth models 34 1. which explains why in that case the inversion problem can be solved, 
since in those models can be negative. 

VI. CALCULATING D L (z) AND mn(z) FOR ACDM MODELS. 

The metric of a ACDM model is the FLRW metric, a special case of LTB solution, where : 

Pair) = Po (60) 
k{r) = (61) 
t b (r) = (62) 
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In this section we will calculate independently the expansion of the luminosity distance and 
the redshift spherical shell mass for the case of a flat AC DM. 

We will also use these formulas to check the results given in the previous section, since in 
absence of cosmological constant and large-scale inhomogeneities they should coincide. 

One of the Einstein equation can be expressed as: 



h(z) = # yOn(^J +n A = H ^/a m (i + z) 3 + n A (63) 

We can then calculate the luminosity distance using the following relation, which is only 
valid assuming flatness: 



"M-v + '^m (64) 



From which we can get: 



Dl = 7T ( 65 ) 

DL = (66) 

It is convenient to re-express the above coefficient in terms of the two observables H , q 

*>l = W (68) 
Dl = ^ (69) 



! _ 3q 2 + g - 2 
' = (i//o 



l ~ —7^ (70) 



where we have used the following relations 

VLl + Qm = 1 (71) 
n M = (72) 

It should be underlined here that H ,q ,ao which appearing in this formulas are not the 
same as the ones defined in the previous section for LTB models. 
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For calculating mn(z) we first need r(z), which can be obtained from the radial null geodesic 
equation which in this case take the simplified form 

dr 1 1 

— = = (73) 

dz (1 + z)d a H 

which can be easily integrated to give 

/ \ 1 f z dz' 

z (q + l)z 2 (3q% + 4q + 1) z 3 



+ lJ(/0 7V J + ( 75 ) 

V\n.r\ H r\ \ ' 



a H 2aoH 6a H 
We can now calculate mn(z) 



p = 3a 6 n M Ht (76) 

/ \ 2 ^4 por(z) 2 dr(z) 

mn{z) = mn,2Z +mn 3 z +mn. = , — - — = (77) 

y/l - k{r{z))r{z) 2 dz 
2(q + l)z 2 4(g + l)V , 5(q + l) 2 (9q + 5)z 



+ ^ + + + O (, 5 ) (78) 



i7 i7 6i? 

We can check the consistency between these formulae and the one derived in the case of 
LTB without cosmological constant by setting: 

k 2 = t\ = (79) 
g = 1/2 (80) 

which corresponds the case in which Qm = 1- 

It is important to mention again that do, qo, H defined in this section are in general different 
from the ones defined in the previous section for LTB models but for simplicity of notation we 
have used the same symbols. 

In the next section, where we will study the inversion problem, we will instead need to clearly 
distinguish between them and for that reason we will introduce Oq , <?o > t° indicate the ones 
corresponding to a ACDM model. 



VII. THE INVERSION PROBLEM CANNOT BE SOLVED 

In this section we will denote with an upper script A all the relevant quantities referred to 
a ACDM model, including the coefficients of the redshift expansion for D L (z) and mn(z). 
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In order to solve he inversion problem we need to solve the following system of six equations 
equations: 

Df = Di l<i<3 (81) 
mnf = mrij 2 < j < 4 (82) 

It is evident that the first and second order term of the luminosity distance expansion imply 
that 

go = g A (83) 
H = < (84) 

If we try match the second and third order term of mn(z) instead we get: 



* = 9o A (85) 

H„ = (86) 
(1 + Qo ) 

It is clear that there is no common solution unless go = \i wn i cn corresponds exactly to set 
to zero the cosmological constant. 

This proof is general, and it does not depend on the particular value of I^a, as long as it is 
non zero. In this regards this is an even stronger statement than the one already well known 
fact that go cannot be negative for a smooth LTB model. 

In fact even if go > 0, which for a flat AC DM would imply Q\ < |, it would still be 
impossible to mimick the cosmological constant for both mn(z) and Dl{z). 

In other words no value of the cosmological constant can be mimicked, even if for positive 
go, which is in principle possible even in a AC DM models with a smaller value of Q\ than the 
presently observed one, which is about 0.7, and which implies a negative go- 

This means that smooth LTB models are fundamentally different from AC DM models, and 
the sign of go it is not the really important observational feature to focus on, but we should 
rather consider the full observational implications at various redshift, since quantities like mn(z) 
encode much more information about the underlying space-time geometry than simply go. 

If we had relaxed the condition t\ — k\ — the inversion problem could be solved, but we 
will give the solution in a separate future work [34( , where we will also try to better understand 
the nature of the singularity occurring in this case. 

From our calculations we have learnt that we can mimick the effects of a cosmological for each 
observable separately, but taking into consideration both of them allow to uniquely distinguish 
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between a smooth LTB and ACDM models. In this sense since an accurate calculation of 
mn(z) from observational data does not exist yet, there is still an open possibility that the 
observed mn(z) is not the one corresponding to the the ACDM model but it is better fitted 
by an appropriate LTB model, in which case mn(z) would be the observable which would allow 
us to rule out ACDM as a realistic cosmological model. 

This was actually the main motivation to introduce Erss{ z )) to test self-consistently in 
redshift space the homogeneity of matter distribution. In this regard Erss{z) is actually more 
suitable for analyzing observational data from galaxy surveys because it can avoid the effects 
of astrophysical evolution of the sources. 

VIII. CONCLUSION 

We have calculated for a central observer in a LTB space without cosmological constant 
the low redshift expansion for the luminosity distance and the redshift spherical shell mass 
respectively to third and fourth order. We have then derived the same observables for ACDM 
models and shown how it is impossible to mimick dark energy with a smooth LTB matter 
dominated universe . Our results imply that future accurate analysis of galaxy surveys data 
could allow us to determine mn(z) to distinguish without any ambiguity between smooth LTB 
and ACDM models. 

What we have shown does not diminish the importance of the study of the cosmological 
implications of large-scale inhomogeneities, but simply give a definite proof that smooth pres- 
sureless spherically symmetric matter inhomogeneities are not sufficient for a central observer 
to explain the observational data corresponding to a flat ACDM model. 
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